Accuracy matrix in generalized simultaneous measurement of a qubit system 
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We formulate the accuracy of a quantum measurement for a qubit (spin-1/2) system in terms of 
a 3 by 3 matrix. This matrix, which we refer to as the accuracy matrix, can be calculated from 
a positive operator- valued measure (POVM) corresponding to the quantum measurement. Based 
on the accuracy matrix, we derive trade-off relations between the measurement accuracy of two or 
three noncommuting observables of a qubit system. These trade-off relations offer a quantitative 
information-theoretic representation of Bohr's principle of complementarity. They can be interpreted 
as the uncertainty relations between measurement errors in simultaneous measurements and also 
as the trade-off relations between the measurement error and back-action of the measurement. A 
no-cloning inequality is derived from the trade-off relations. Furthermore, our formulation and 
the results obtained can be applied to analyze quantum-state tomography. We also show that the 
accuracy matrix is closely related to the maximum-likelihood estimation and the Fisher information 
matrix for a finite number of samples; the accuracy matrix tells us how accurately we can estimate 
the probability distributions of observables of an unknown state by a finite number of quantum 
measurements. 

PACS numbers: 03.67.-a, 03.65.Ta, 03.65.Wj 



I. INTRODUCTION 

Accessible information about a quantum system is re- 
stricted by the noncommutability of observables. The 
nature of this restriction can be classified essentially into 
two categories: fluctuations inherent in a quantum sys- 
tem and the error caused by the process of measurement. 
These aspects of uncertainty constitute the two distinc- 
tive features of quantum mechanics. 

The Kennard-Robertson uncertainty relation such as 
AzAp > h/2 describes quantum fluctuations that are 
independent of the measurement process [E H H HI- 
According to Bell's theorem Q, this type of quantum 
fluctuations prohibits us from presupposing any "ele- 
ment of reality" [6j behind the probability distributions 
of observables. The measurement error, on the other 
hand, is determined by the process of measurement which 
is characterized by a positive operator-valued measure 
(POVM) 0,0. In the idealized error- free limit, quantum 
measurement is described by projection operators which, 
however, cannot always be implemented experimentally 

Information about more than one observable can be ob- 
tained from a single POVM in simultaneous measurement 
of two noncommuting observables and quantum-state to- 
mography. It is known, however, that, in simultaneous 
measurements, at least one of the observables cannot be 
measured without incurring a measurement error [9(. In 
this context, various uncertainty relations between the 
measurement errors of noncommu ting observables have 

been studied [H El El El El EE M ES El El H HI 

In this paper, we quantify the measurement accuracy 
and the measurement error of observables in terms of 
a given POVM E = {Ek} by introducing 3x3 accu- 



racy matrix x(E) calculated from the POVM. Based on 
this accuracy matrix, we derive trade-off relations be- 
tween the measurement accuracy of two or three observ- 
ables, these being stronger trade-off relations than those 
derived in our previous work [23|. They can be inter- 
preted as the uncertainty relations between the measure- 
ment errors of noncommuting observables in simultane- 
ous measurements or as the uncertainty relations between 
the measurement error and back-action of the measure- 
In addition, a no-cloning in- 
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equality 3(J" 3ll 32L |33| is derived from the trade-off re- 
lations. 

In a rather different context, the maximum-likelihood 
estimation [34], [HJ has been investigated as the standard 
scheme of quantum state tomography for a finite num- 
ber of samples. Several studies have focused on the ef- 
ficiency and optimality of the estimation of an unknown 
quantum state [H, [33, [H, [U HJ. We show that our 
characterization of the measurement accuracy can be re- 
lated to the maximum-likelihood estimation and that the 
accuracy matrix can be interpreted as an average of the 
Fisher information matrix over the state to be measured. 
The trade-off relations can also be interpreted as those 
concerning the accuracy of the estimate of various prob- 
ability distributions of noncommuting observables. 

The constitution of this paper is as follows. In Sec. I 
I, we formulate the general quantum measurement of a 
qubit (spin-1/2) system. In Sec. Ill, we define the ac- 
curacy matrix and investigate its properties. Based on 
this accuracy matrix, we define the accuracy parameter 
and error parameter in a particular direction of measure- 
ment. In Sec. IV, we derive the trade-off relations be- 
tween the accuracy parameters or the error parameters in 
two or three directions. In Sec. V, we apply the trade-off 
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relations to specific problems: the uncertainty relations 
between measurement errors in nonideal joint measure- 
ments, the uncertainty relations between the error and 
back-action, a no-cloning inequality, and quantum state 
tomography. In Sec. VI, we point out a close connection 
between the accuracy matrix and the Fisher information 
matrix. We conclude this paper in Sec. VH. 

II. QUANTUM MEASUREMENT OF A QUBIT 
SYSTEM 

We consider a quantum measurement described by 
POVM E = {E k } (k = 1, 2, • • • , m) on state p of a qubit 
system, where k denotes the outcome of the measure- 
ment. POVM E satisfies YJ fc E k — I, with / being the 
identity operator, and can be parameterized as 

E k = r k (I + v k -&), (1) 

where & = (<J x ,<7 y ,cr z ) represents the Pauli matrices. 
The requirements that the sum of E k s equals the iden- 
tity operator and that all of them be nonnegative are met 
if and only if 

J^r k = 1, rkVk = °' Tk > °' \ Vk \ - 1 for a11 k - ( 2 ) 

fc k 

We can also parameterize density operator p as 

p =l(i + s -a), (3) 

where Sq is the Bloch vector satisfying |so| < 1. Con- 
versely, for a given />, Sq is calculated as Sq — tr(pcr). 
The probability of obtaining the measurement outcome 
k is then given by 

q k = tr(E k p) = r k (l + v k ■ s ). (4) 

Any observable O of the qubit system can be diago- 
nalized as 

6 = \+P(+;n) + A_P(-;n), (5) 

where A+ and A_ are the corresponding eigenvalues, 
P(+;n) and P(+;n) are projection operators with n 
being a three-dimensional unit vector, and 

P(±;») = i(/±»-*). (6) 

The probability distribution of observable O is then given 
by 

p(±; n) = tr(P(±; n)p) = ~(1 ± n • s Q ). (7) 

If we are not interested in eigenvalues of the observ- 
ables but are only concerned with the directions (±) of 
the outcome, we can replace O with n ■ & by setting 



\± = ±1. In the following analysis, we identify observ- 
able A + P(+;n) + A_P(— \n) with the observable n ■ & 
and refer to the probability distribution in Eq. ([7]) as that 
in the direction of n. 

We discuss three typical examples. 

Example 1 (projection measurement). We can pre- 
cisely measure n • <r by the projection measurement de- 
scribed by the POVM E = {P(+;n),P(-;n)}. 

Example 2 (Nonideal measurement). A more general 
class of measurements can be described by the POVM E 
consisting of two positive operators parametrized as 

E(+\n) = r(I + e 1 n-&), E(-;n) = (1 - r)(J- e 2 n- &), 

(8) 

where n is a unit vector, re\ — (1 — r)e 2 = 0, < r < 1, 
-1 < e x < 1, and -1 < e 2 < 1- This POVM cor- 
responds to a nonideal measurement of the observable 
n ■ & [H, It can be reduced to a projection measure- 
ment {P(+; n), P(— ; n)} if and only if e± = e% = 1 and 
r = 1/2. On the other hand, the POVM is trivial (i.e., 
E + = rl and E- = (1 — r)I) if and only if S\ = £2 = 0; 
then we cannot obtain any information about p. Equa- 
tions © can be rewritten as 

where F is the 2x2 transition-probability matrix 

T?-( r (! r(l-ei) \ nm 

"l,(l-r)(l- E2 ) (l-r)(l+e 2 ) J 

which satisfies J2i F ij = 1 and < |detF| 2 < 1. Note 
that F describes a binary symmetric channel ^4lj if and 
only if r = 1/2 and e\ = e 2 - It follows from Eq. © that 
any measurement process described by a POVM consist- 
ing of two positive operators is formally equivalent to a 
measurement process in which a classical error is added 
to the projection measurement. The physical origin of 
this error, however, lies in the quantum-mechanical in- 
teraction. 

Example 3 (probabilistic measurement). Suppose 
that a nonideal measurement of A = tia ■ <J is performed 
with probability £ (0 < £ < 1) and that B = tib ■ a 
is performed with probability 1 — £. The POVM cor- 
responding to this probabilistic measurement consists of 
four operators: 

E = {£E(±;n A ),(l-Z)E(±;n B )}. (11) 

As the number of measured samples increases, this mea- 
surement asymptotically approaches the measurements 
on ./V identically prepared samples which are divided into 
two groups in the ratio 6 : 1 — £, with A being measured 
for the first group and B for the second group. 

Other important examples such as nonideal joint mea- 
surements and quantum state tomography are discussed 
in Sec. VI. 
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III. ACCURACY MATRIX 
A. Definition of the Accuracy Matrix 

We will characterize the accuracy of an arbitrary ob- 
servable in such a manner that it depends only on the 
process of measurement and not on the measured state 
p. We first define the accuracy matrix. 

Definition 1 (accuracy matrix). The 3x3 accuracy 
matrix x(E) characterizing the measurement accuracy of 
observables in terms of the POVM E is defined as 

x( E )ij = J2r k (v k ) l (v k ) j , (12) 

k 

where (vk)i denotes the ith component of the real vector 
v k and ij shows indices of matrix elements of x(E). We 
introduce the notation vv T with v £ K 3 as 

(vv T ). l] = (vUv) r , (13) 

that is, v T denotes the transposed vector of v and vv 1 
denotes the projection matrix onto direction v in K 3 
whose ij matrix element is given by (v)i(v)j. We can 
then rewrite (|12|) in matrix form as 

x(^) = J2 TkVkv k- ( 14 ) 

k 

Note that x( E ) is positive semidefinite and Hermitian, 
and can therefore be diagonalized by an orthonormal 
transformation. 

The physical meaning and useful properties of the ac- 
curacy matrix will be investigated subsequently, and its 
foundation from an information-theoretic point of view 
will be established in terms of the maximum-likelihood 
estimation of the probability distribution of observables 
in Sec. VI. In fact, the accuracy matrix is closely related 
to Fisher information matrix p07p or (|108[) . although 
physical quantities such as the measurement error can 
be directly derived from the accuracy matrix without re- 
sort to Fisher information. 

Noting that J2k rk \ Vk \ 2 — ^lk Vk = 1> we can obtain 
the following fundamental inequality which forms the ba- 
sis of trade-off relations to be discussed later. 

Theorem 1. Three eigenvalues {xi,X2,X3j OI x(E) 
satisfy 

Xi+X2+X3<l, (15) 

or equivalently, 

Sp( X (E)) < 1, (16) 

where we denote the trace of the 3x3 matrix as Sp(- • • ) 
to reserve symbol tr(- • • ) for the trace of a quantum- 
mechanical 2x2 matrix. The equality xi + X2 + X3 — F 
or Sp(x) = 1, holds if and only if \v k \ = 1 for all k. 

The following corollary follows from the positivity of 
X(E). 



Corollary 1. The accuracy matrix satisfies the fol- 
lowing matrix inequality: 

< X(E) < I 3 , (17) 

where -Z3 is the 3x3 identity matrix, and x( E ) < I3 
means that all eigenvalues of I3 — x( E ) are nonnegative. 

The following examples illustrate the physical meaning 
of the accuracy matrix. 

We first consider a nonideal quantum measurement 
(see also example 2 in Sec. II). We can rewrite Eq.© 
as 

Eh. = r(I + V! ■ &), E 2 = (l- r)(J + v 2 ■ &), (18) 

where V\ = e\n and v 2 = —£2^- The accuracy matrix 
can then be represented by 

X(E) ee rvivj + (1 - r)v 2 vj = xn™n T , (19) 

where xn is th e eigenvalue of x corresponding to the 
eigenvector n, and is given by 

X11 = r\ Vl \ 2 + (1 ~ r)\v 2 \ 2 =eie 2 . (20) 

We can also write xn m terms of the transition- 
probability matrix introduced in Eq. (|10p as 

_ |detJf IdetJf _ Ideti^l 2 
XU 4r 4(1 -r) 4r(l - r) ' 1 ' 

The accuracy parameter xn satisfies 

< X11 < 1, (22) 

where xn — 1 holds if and only if \vi\ = \v 2 \ = 1 and 
r = 1/2; that is, E describes the projection measurement 
of observable n ■ <x. Note that x(E) = nn T holds in 
this case. On the other hand, xn = holds if and only 
if |vi| = \v 2 \ = 0. In this case, x(E) = O holds, and 
we cannot obtain any information about p. The nonzero 
eigenvalue xn thus characterizes the measurement ac- 
curacy of n ■ <t; the larger xn> the more information 
we can extract about n ■ <x from the measurement out- 
come. These properties can be generalized for an arbi- 
trary POVM as shown below. 

Another example is the probabilistic measurement of 
two noncommuting observables (see example 3 in Sec. I 
I). We consider the nonideal measurement of A whose 
accuracy matrix is XATiAn^ and that of B whose accu- 
racy matrix is XB n B n ]j- The accuracy matrix of the 
probabilistic measurement is given by 

X(E) = ZxAUAnl + (1 - OxbUbuI. (23) 

This representation suggests that the measurement accu- 
racy concerning A is degraded by a factor of £ compared 
with the single nonideal measurement of A, because we 
cannot observe A with probability 1 — £. A similar ar- 
gument applies to B as well. Equation (f2"3")l shows that 
x(E) is the linear combination of the accuracy matrices 
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of POVMs measuring A and B, where the coefficients £ 
and f — £ give the probabilities of measuring A and B, 
respectively. 

This can be generalized as follows. Let us con- 
sider three POVMs: E' = {E u E 2 , ■ ■ ■ , E m }, E" = 
{E m+ i, E m+ 2 7 • • • , E n }, and E = {£Ei, • • • , £E m , (1 — 
£)E m+ i, ■ • • , (1 - E n} with < £ < f . The POVM E 
describes the probabilistic measurement of E' with prob- 
ability £ and that of E" with probability 1 — £. According 
to the definition of the accuracy matrix, we have 



X (E) = Y,(Zr k )v k vJ + {(1 - 0»*}«k«fc 

fc=l k—m-\-l 

= £t(E') + (1 - Ox(E"). 

We thus obtain the following theorem. 
Theorem 2 (linearity): 

x(E)=€x(E') + (l-Ox(E"), 
or more symbolically, 

X (£E' + (1 - OE") = Sx(E') + (1 - C)x(E"). 



(24) 



(25) 



(26) 



Note that we can take as a scalar measure of the mea- 
surement accuracy the largest eigenvalue of the accu- 
racy matrix which we denote as x(E) max . It satisfies 
< x( E )max < 1, where x(E) max = 1 holds if and only if 
E describes the projection measurement of a particular 
direction and x(E) m ax = if and only if the POVM 
is trivial: E = {q k I}, where / is the identity opera- 
tor and q k denotes the probability of finding outcome 
k, with ^2 k qk — 1- We may alternatively choose the 
scalar measure to be Sp(x(E)); it has the linear property 
from theorem 2 and satisfies < Sp(x(E)) < 1, where 
Sp(x(E)) = if and only if the POVM is trivial. 



B. Accuracy Parameter in a Specific Direction 

We next parametrize the measurement accuracy of a 
particular observable. We denote the support of x(E) as 
V(E); that is, V(E) is the subspace of R 3 spanned by all 
eigenvectors of x(E) with nonzero eigenvalues. 

Definition 2 (measurement accuracy). The accuracy 
parameter x( n ; E) in direction n <E V(E) is defined as 



x(n;E) 



n.(x(E)->' 



(27) 



where x(E) 1 is assumed to act only on subspace V(E). 
If n ^V(E), we set x(n; E) = 0. 

This definition is closely related to the Fisher in- 
formation concerning a particular direction defined in 
Eq. (fTO . 

Definition 3 (measurement error) The error param- 
eter of the measurement in direction n is defined as 



e(n;E) 



1 



x(n;E) 



l = n-(x(E)- 1 )n-l. 



(28) 



The parameters x( n l E) and e(n; E) satisfy the follow- 
ing inequalities. 
Theorem 3: 



o<x(«;E) < l, 

< e(n;E) < oo. 



(29) 
(30) 



The equality x( n ; E) = 1, or cquivalently e(n;E) = 0, 
holds if and only if the measurement described by E is 
equivalent to a projection measurement in direction n. 
In this case, the measurement involves no measurement 
error. The other limit of x( n i E) = 0, or equivalently 
e(n; E) = oo, holds if and only if n £V(E). In this case, 
we cannot obtain any information about direction n from 
the measurement. 

Proof Since x(E) commutes with the identity oper- 
ator ^3, we can show that I3 < x(E)" 1 from inequality 
(fT7| in corollary 1. We thus obtain 



1 = n ■ I 3 n < n ■ ( X (E)' 



(31) 



Inequalities (|29|) and (|30|l are the direct consequences 
of this inequality. The condition that x( n ! E) = and 
e(n; E) = 00 hold follows from the definitions of x( n ; E) 
and e(n; E). 

We next show the condition that x( n ! E) = 1 and 
e(n; E) = hold. If E is the projection measurement 
in direction n, then x( ri ;E) = 1. Conversely, from in- 
equality (fT5|) . it can be shown that if x^jE) = 1 and 
e(n; E) = hold, then n is the eigenvector corresponding 
to eigenvalue 1 and that the other two eigenvalues are 0. 
It follows from the condition of equality Sp(x(E)) = 1 in 
theorem 1 that v k = 1 for all k. Therefore, without loss 
of generality, we can write the POVM as 

Ex =n(I + n ■&),■■■ , E m =r m (I + n-a), (32) 
E m+1 = r m+ i(I -n- &),■■■ , E n = r n (I -n-&), (33) 

where r k = £Lm+i r k = 1/2, because £ fe r k = 1 

and Y^kLi r k n — J2k= m +i r k n = hold. We define two 
operators as 



P(+;n) 



fc=i 



E k , P(-;n) 



E 

k—m-\-l 



E k : 



(34) 



then {P(+; n), P(— ; n)} describes the projection mea- 
surement in direction n. — 

Let Tii, ri2, and 713 be the eigenvectors of x(E), and 
Xi, X2, and \3 be the corresponding eigenvalues. It can 
be shown that 



x(n i; E) = x*, (i = l,2,3). 



(35) 



According to theorem 1, we cannot simultaneously mea- 
sure the three directions corresponding to the eigenvec- 
tors with the maximum accuracy x» = 1 for all i; the 
trade-off relation (fT~5|) or (flB")) is equivalent to 



Xi +X2 + X3 < 1. 



(36) 
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This trade-off relation represents the uncertainty relation 
between the measurement errors in the three directions. 

We define that the POVM E is optimal if and only if 
Sp(x(E)) = xi + X2 + X3 = 1; that is, E reaches the 
upper bound of trade-off relation (fl~5|) , (fl6|) . or (f36|) . On 
the other hand, we define that E is symmetric if and only 
if E) = x{ n '\ E) holds for any n and n' . In this case, 
x(E) is proportional to the 3x3 identity matrix. 



C. Reconstructive subspace 

We next introduce the concept of "reconstructive sub- 
space" and "reconstructive direction." The following the- 
orem can be directly shown from the definition of the 
accuracy matrix. 

Theorem 4 V^(E) corresponds to the subspace 
spanned by the set of basis vectors {v k } of the accuracy 
matrix (fT4| . 

Suppose that we perform the measurement {E k } and 
obtain the probability distribution {qk} for each outcome 
k. Can we then reconstruct the premeasurement distri- 
bution {p(n;±)} of the system from {(ft}? The answer 
is given by the following theorem. 

Theorem 5 (reconstructive subspace and reconstruc- 
tive direction) We can reconstruct the probability distri- 
bution {p(±;n)} from the measured distribution {q k } if 
and only if n 6 V^(E). We thus refer to V(E) as a re- 
constructive subspace and to a unit vector in T^(E) as a 
reconstructive direction. 

Proof We can show from Eq. ((4]) that 




(37) 



where M is a m x 3 matrix 



V w 



M 



( n(vi) x n(vi)y 

r 2 (v 2 ) x r 2 (v 2 ) y 



ri(vi) 2 
r 2 (v 2 ) z 



(38) 



y Tm{Vm)x ^ra('^m)y ^mi^ra)z J 



Let Ker(M) and Im(M) be the kernel and image of M, 
respectively. It can easily be shown that Ker(M) = 
^(E)^. Let us introduce the equivalence relation "~" 
as a — 6 <^ a - b e V^E)^. We denote the 
equivalence class of v S R 3 as [v], where [v] is an el- 
ement of the quotient space R 3 / ~. From the homo- 
morphism theorem, the quotient map M/ ~ is a lin- 
ear isomorphism from R 3 / ~ to Im(M). Noting that 
(<Zi,92, • ■ -qm)* - {n,r 2 , ■ ■ -r m )* € Im(M), we obtain 



[ So ] = (M/ 



/ 








\ 




92 




r 2 




V 






\r m ) 


1 



(39) 



By taking a representative s' £ [so], we can reconstruct 
n ■ s as n ■ s = n ■ s' which gives p(±;n) through 
Eq. Q.- 

We consider the nonideal measurement of A — n A ■ & 
with the POVM E A = {E(+; n A ), E(-;n A )} in Eq. ©. 
The nonideal measurement is characterized with the ac- 
curacy matrix x^n^n^. In this case, we can show that 
V(E A ) = {a | a = \n A ,\ e R}. It follows that 
x(±n A ;~E) = xa and x( n i E) = for n ^ ±n A . 

We next consider the probabilistic measurement of E^ 
and Eg in example 3 in Sec. H. The probabilistic mea- 
surement is characterized with the accuracy matrix of 
the joint POVM E given in (|23|) . so the reconstruc- 
tive subspace is two dimensional: V(F,) = {\ A n A + 
^B n B I (^a,^b) G R 2 }. A straightforward calculation 
shows that 



x( n A', E) = 064, x(nB] E) = (1 - £)xs- 



D. Data Processing Inequality 



(40) 



If the classical noise described by a transition- 
probability matrix is added to the measurement out- 
comes, the measurement accuracy should deteriorate. 
This fact can be expressed as a data processing inequal- 
ity. 

Theorem 6 (data processing inequality). Suppose 
that two POVMs E = {E x , E 2 , ■ ■ ■ , E n } and E' = 
{E[, E 2 , ■ ■ ■ , E' m } are related to each other by 



(41) 



k=i 



where is an m x n transition-probability matrix sat- 
isfying Y^j Fjk = 1. It then follows that 



X(E') < X (E), 



(42) 



where matrix inequality (|42[) means that all the eigenval- 
ues of x(E) — x(E') are non-negative. 

Proof We can parametrize the POVMs as 



E k =r k (I + v k - a), 



<7 , 



(43) 



(44) 



where = Fj k r k . Introducing the function /„ 
R 3 , with arbitrary vector w G R 3 as 

f w (v) = w ■ (vv T )w = \v ■ w\ 2 , 

we can show that 



w ■ x(E)w = X! r kfw(v k ) 



(45) 



(46) 



fc=i 
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,•=1 \fc=l / v ^fe=l 



(47) 



The Hessian of f w , which is defined as H(f w )ij = 
d 2 f w ( v ) /d{v)id(v)j , becomes 

H(f w ) = ww T > 0, (48) 

so that f w is a concave function. Therefore 



fw i^tkVk j < 

V fc / k 



K) 



(49) 



holds for any satisfying ^ fe = 1 and < £fc < 1 
for all fc. Taking £ fc = Cjfe/Sfe'=i ^'fe') inequality (05]) 
becomes 



(»*)• (51) 



or equivalently, 



\k=l ) V 2^=1 J k=1 

Noting that X^'=i = rfe ' we obtain 

j=l \fc=l / V 2^k=l^k J k=1 



which implies that 

w ■ x(E*')w < w ■ x(E)tt). 



(52) 
(53) 



Since (|53")) holds for arbitrary w, we obtain (|4"2"|) . — 

The following corollary is a direct consequence of the 
foregoing theorem. 

Corollary 6 Suppose that E' = {E[,E 2 ,--- , E' m } is 
obtained by a coarse graining of E = {E\, E 2 , • • ■ ,E n }: 
E[ =E 1 +E 2 + ---+E i(1) , E' 2 = E m+1 + ---+E im , 
and E' m = E i(m _ 1)+1 + ■ ■ ■ + E n , with 1 < < i(2) < 
■ ■ ■ < i(m — 1) < n. Then 



X(E') < X (E) 



(54) 



holds. Inequality ([54")) means that the measurement ac- 
curacy in any direction is decreased by a coarse graining. 

We can also express the data processing inequality in 
terms of the accuracy parameter in an arbitrary direc- 
tion. 

Theorem 7 (data processing inequality) . We consider 
the POVMs E and E' satisfying Eq. (gl]). Suppose that 
V(E) = V(E') = R 3 . Then, 



^(^-^^(E)- 1 
holds, or equivalently, 

x(n;E')<xO;E) 



(55) 
(56) 



holds for arbitrary n. 

Proof Let xij X2, and X3 be the eigenvalues of x(E), 
and rti, n 2 , and rj.3 be the corresponding eigenvectors. 
Similarly, let x'i, X21 and X3 be the eigenvalues of x(E'), 
and n[, n' 2 , and n' 3 be the corresponding eigenvectors. 
It follows from the data processing inequality (|4"2")l that 

3 

«. • x(E')ni = n irfj ^Xi = rii- x(E)n,, (57) 

i=i 

for i = 1, 2, 3, where Jiy = rii ■ n'j. Applying the concave 
inequality to I fx, we obtain 



x- x < E^;- • »<w '• ( 58 ) 

For arbitrary n, we can show that 

3 



n-(x(E)- 1 )n = E( n ' n ') 2 ^ 

i=l 
3 



2^2 „,/-i 
3 



i,3=l 

jyi-^x 1 - 1 

n-( X (E')- 1 )n, 



(59) 



which implies (|55|) and ([56")) . 



IV. TRADE-OFF RELATIONS FOR 
GENERALIZED SIMULTANEOUS 
MEASUREMENT OF A QUBIT SYSTEM 

We now derive general trade-off relations between the 
measurement errors of noncommuting observables, which 
are the main results of this paper. 

Let ni, n 2 , and n.3 be the respective eigenvectors of 
x(E) corresponding to the eigenvalues xij X2, and X3i 
where \i = x( n i;E) (i = 1,2,3). We define the error 
parameters as £; = £(n^;E) = (xi) — 1 — 1- Inequality 
(|15p or (|16p in theorem 1 can be rewritten in terms of 
the error parameters as 



£i£2£3 > £1 + £2 + £3 + 2. 



(60) 



Considering two eigenvalues alone (i.e., Xi + X2 < 1), we 
can simplify the trade-off relation: 



£i£2 > L 



(61) 



The trade-off relations (fBT)]) and ([BT]) can be generalized 
to the case of arbitrary directions. We first consider the 
case of two observables. 

Theorem 8 (trade-off relation). We consider a si- 
multaneous measurement in two directions i%a and ng 
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(tia -tib = cos 9) described by the POVM E. We assume 
tia € V(E) and n B e V(E), and define e a = e(n a ;E) 
and Xa = x( n a', E) (a = A,B). Then the trade-off rela- 
tion 

e A e B >sm 2 9, (62) 

or equivalently, 

Xa + Xb ~ XaXb cos 2 6 < 1 (63) 

holds. 

Proof. We divide the proof into two steps. 

Step 1. We consider a situation in which both tia and 
n B lie in a plane spanned by two eigenvectors. With- 
out loss of generality, we choose ri\ and ti 2 as the two 
eigenvectors, and expand ua and % as 

tia — ni cos 9a + Ti2 sin 9a, (64) 

n.B = ni cos#b + n.2 sinflg, (65) 
where 9 — 9 'a — 9b- It can be shown that 

ea = £1 cos 2 9a + £2 sin 2 9a, (66) 



Combining this with J^i Qij = 1, we obtain 

3 n 3 3 n 2 3 i 

I— 1 6 2 — 1 J — 1 J J — 1 17 

Therefore 

1 + 1 <V— ^- <V 1 <!■ (74) 

1 * l—l 1 J—l J 

This inequality means that x'i + X2 — 1) or equivalently, 
44 > 1. (75) 

We can derive inequality (|62|) by following the same pro- 
cedure as in step 1. We can directly derive inequality 
p?)) lrom (EH)-— 

We note that the equalities in (|6"2"|) and (|63[) hold in 
the case that the POVM E is given by {|aij| j ± cc, • <r} 
(i = 1, 2), where Xi = (^xati^ + ^xsn B )/4 and x 2 = 
(Vxa^a - y/xB~n B )/4:. 

The accessible regime for \A and Xb is illustrated in 
Fig.l for the case of d = vr/2, 9 = tt/6, and 9 = 0. Note 
that regime Q can be reached only th roug h simultaneous 
measurement for the case of 9 = ir/6 [231 ] - 



cb = £i cos 9 B + £2 sin 



(67) 



Applying the Cauchy-Schwarz inequality and £162 > 1, 
we obtain 



£a£b > (y/e~i£2(cos a sin 9 B 
= £i £2 sin 2 
■> sin 2 9. 



sm I 



icos0 s )) 2 (68) 
(69) 
(70) 



The equality sa£b — sin 9 holds if and only if 
E\ cos 9a cos #b + £2 sinf?^ sin 6>s = and £i£2 = 1. 
In the case of sa = £s (i- e -j the measurement er- 
rors are symmetric), the equality holds if and only if 
sin(6U + B ) sin(9 A - 9b) cos(9 a + 9 B ) = 0. 

Step 2. We next consider a more general case. We 
choose an orthonormal basis {n[, n 2 ,n' 3 } such that both 
tia and n B are in the plane spanned by n' x and n' 2 . 
We introduce the notation = e(ti^;E) (i = 1,2,3). 
Let Qij be a 3 x 3 orthogonal matrix which transforms 
{ni, TI2, 113} into {ti'i, n' 2 , n' 3 }. It can be shown that 



2 -£- 



(71) 



Note that £\ Qf^ = 1 because Qij is an orthogonal ma- 
trix, and that the function (1+x) 
from a concave inequality that 



is concave. It follows 



1 h-> 



(72) 



0.5 



0.6 



rr. 
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FIG. 1: Trade-off relation for the accuracy of noncommuting 
observables. P indicates the regimes satisfying the inequality 
for the case of 6 — it/2, the union of P and Q indicate the 
regime satisfying inequality (|63p for the case of 9 = tt/6, and 
the union of P, Q and R indicate the regime satisfying the 
inequality for the case of 9 = 0. We can only access regime Q 
through simultaneous measurement for the case of = tt/6. 

The trade-off relation can be interpreted as the uncer- 
tainty relation between measurement errors [9j, MM, 

MMMMMMMMMMMMm- it offers 

a rigorous representation of Bohr's principle of comple- 
mentarity [l0( which dictates "the mutual exclusion of 



8 



any two experimental procedures" when we measure two 
noncommuting observables simultaneously. 

The trade-off relation between three observables can 
be formulated as follows. 

Theorem 9 We consider a simultaneous measure- 
ment in three directions n A , fig, and nc described by 
the POVM E. Let us assume that n A , n B , and nc are 
linearly independent. We set the notation e a = £{n a ; E) 
and Xa = x( n ctiE), where a = A,B,C. Then the in- 
equality 



£a£b£c > %{n A • (n B x n c )} 2 



(76) 



holds. The equality in ([76]) holds if and only if E\ — £2 — 
£3 = 2 and {n a } are orthogonal. 
Proof Introducing the notation 



n„ 



£2(n a ) 2 
£3(n a ) 3 



(77) 



where (n a )i = n a ■ ra, (a — A,B,C, i — 1,2,3), it can 
be shown that 



(78) 



We thus obtain 



£a£b£c = \h A \ 2 \h B \ 2 \h c \ 2 



> {h A ■ (n B x h c )} 2 

= £i£ 2 £z{n A ■ (n B x n c )} 2 

> (£1 + £2 + £3 + 2){n A ■ (n B x n c )} 2 . 



(79) 



We can show e\ + £2 + £3 > 6 from Xi + X2 + X3 < lj 
therefore we obtain ([751) • — 



V. APPLICATIONS 

We have discussed in Sec. IV trade-off relations (|f)2"]) 
and (|63[) which describe the uncertainty relations in gen- 
eralized simultaneous measurements. In this section, we 
discuss possible applications of these trade-off relations. 



A. Nonideal Joint Measurement 

We consider a class of simultaneous measurements 
called nonideal joint measurements, where two observ- 
ables A = i%a ■ (T and B = n B ■ & are simultaneously 
measured. Since their eigenvalues are ±1, each measure- 
ment should give a pair of outcomes (i,j) — ±) for 
observables A and B. The joint POVM E = {E(i,j)} 
can be parametrized as 



E{i,j) = rij(I + Vij ■ <x). 



(80) 



The marginal POVMs E Q = {E a {i)} {a = A, B) are 
defined by 

E A {+) = E(+, +) + E(+, -), E A {-) = 1- E A (+), 

E B (+) = E(+,+) + E{-,+), E B (-) =1- E B {+), 

(81) 



and can be parametrized by 

E a (+) = r a (I + v a ■ <t), 

where 

r++v+++ r+-v+- 



r A = r ++ +r+_ , v A 



r B = r ++ + r_ + , v B = 



; - +V++ + r_ + t;_ 
r++ +r~+ 



(82) 



(83) 



This simultaneous measurement can be regarded as a 
nonideal joint measurement (lol . Il7| if and only if the 
marginal POVM E A ( B j corresponds to the nonideal mea- 
surement of A(B), that is, 

v a \\n a . (84) 

We can define the 2x2 transition-probability matrices 
of E Q as in Eq. (TTU]): 



F a = 



r a (l±\v a \) r a (l =F \v a \) 
1 -r a (l± \v a \) l~r a (lT\v a \ 



(85) 



In this case, we can calculate the accuracy of A and B 
by two different methods. One method to calculate the 
accuracy parameter is based on the joint POVM E: 



x Q = x(™ Q ;E), 



(86) 



where a = A,B. The other is based on the marginal 
POVM E Q : 



xL = x(™ Q ;E Q ). 



(87) 



Note that x( e q) = x' a n an^. 

These two accuracy parameters are equivalent as 
shown by the following theorem. 

Theorem 10: 



Xa Xo 



(88) 



The proof of theorem 10 is given in the Appendix. 

For n 7^ ±n a , we can show that n is not an element 
of V(E a ); therefore, x( n ;E Q ) = 0. On the other hand, 
x( n ', E) > holds by definition. We can thus obtain the 
following corollary. 

Corollary 10 For arbitrary n, 



x(n;E Q ) <x(n;E). 



(89) 



We next discuss the relationship between the present 
work and our earlier work 12311 for the case of nonideal 



9 



joint measurement. In Ref. [231 ], we have introduced the 
accuracy parameter X a and error parameter £ a as 

X a = (detF a ) 2 , £ a =l--i. (90) 

On the other hand, the accuracy parameter x' a and error 
parameter e a in the present paper are given by 



Xo 



Xa 



(detF Q ) 5 



Mi -*•«)' 

It can be easily shown that 

Xa __ *^-a i £"<_ !_: j 



l 

1. 



(91) 



(92) 



so the trade-off relations derived in the present paper 
are stronger than our previous ones (£a£b _ sin 2 9 and 
+ A s - X A X B cos 2 < 1) derived in Ref. (23] . The 
latter trade-off relations can thus be derived from those 
obtained in the present paper. 



Uncertainty Relation between Measurement 
Error and Back-action 



We have interpreted trade-off relation (I62[) as the un- 
certainty relation between the measurement errors. In 
this subsection, we show that it can be interpreted as the 
uncertainty relation between the measurement error and 
back-action of the measurement [H [H, [U _S __] . Let 
us suppose that p' is a state immediately after the mea- 
surement of A = tia ■ <r for the premeasurement state p. 
If the measurement of A is described by measurement op- 
erators {Mfc}, we can write p' as p' = Ylk^kfiM'^ ■ F° r 
simplicity, we assume that the number of measurement 
outcomes is 2: k — 1, 2. 

To identify the disturbance of B = ng • a caused 
by the measurement of A, we consider how much in- 
formation about B for the premeasurement state p re- 
mains in post-measurement state p'. We character- 
ize this by considering how much information on p 
can be obtained by performing the projection mea- 
surement of B for p' . Note that we can regard the 
projection measurement of B on p' described by the 
POVM as the measurement of p described by the POVM 
{Ek^kM+Wk, EkMiP B (-)Mk}, where n B • <x = 
Pb(+) — Pb(—)- The joint operation of measurement A 
followed by measurement B can be described by a POVM 
{E{i,j)}, where 

E(i,j) = MtP B (j)Mi. (93) 
We can construct the marginal POVMs as 

^E(i,j) = MjMi, Y,E(i,j) =Y,MlP B (j)M k . 

j i i 

(94) 

It is possible to interpret 1 — xb as a measure of the 
back-action of B caused by measurement of A. Defining 



the measurement error of A as ea = (1/xa) — 1 an d the 
back-action of the measurement on B as <Ib = (1/Xb) - L 
we can obtain the trade-off relation between the error and 
back-action based on inequality (f()2"|) . 

Theorem 11 (uncertainty relation between measure- 
ment error and back-action). 



£acLb > sin 6 



(95) 



We note that a non-selective measurement process for 
A can simulate the decoherence caused by the environ- 
ment. In this case, the trade-off relation (pj)) gives a 
lower bound on the back-action of B in the presence of 
decoherence characterized by xa- 



C. No-cloning Inequality 

Another application of the trade-off relation is the 
derivation of a no-cloning inequality. We consider a 
quantum cloning process from qubit system P to qubit 
system Q described as follows: Let p be an unknown 
density operator of system P to be cloned, po be that 
of system Q as a blank reference state, and p cnv be 
that of the environment. The density operator of the 
total system is initially given by p ® po (8 Penvt an d 
becomes Up £g> po (g> p en vU' after unitary evolution U. 
We define p P = trQ :Cnv (Up ® p <g> /w^) and pQ = 
trp,cnv(Up (8 po (E> PcnvU^). We can write pp and pq 
in the operator-sum representation as pp = J2k M k pMl 
and P Q = E k M k pM'l _____ 

The no-cloning theorem [3_, _l|, ]3j^ states that there 
exists no unitary operator U that satisfies pp = pq = p 
for arbitrary input state p. If U is the identity oper- 
ator, then all information about p remains in system 
P, and no information is transferred into system Q; 
pp = p and pq = po- As another special case, if U 
describes the swapping operation between P and Q (i.e., 
Up®po® pcnvU^ — po <8> p ® p cnv ), then all informa- 
tion about p is transferred into Q with no information 
left in P. Intermediate cases between the identity oper- 
ation and the swapping operation can be quantitatively 
analyzed by the no-cloning inequality [33|. 

We derive here another simple no-cloning inequal- 
ity based on the trade-off relation. We first consider 
how much information about p remains in pp. We can 
characterize this by considering how much information 
about n ■ & of p can be obtained by the measurement 
of n ■ it on pp. We can regard the measurement of 
n ■ a on pp as the measurement described by the POVM 

Ep(t») - {E fe &lH+\ n)M k , E fe m£p(-; n)M k } on p, 
where n & = P(+; n) — P(— ; n). We can thus character- 
ize the amount of information that remains in P by the 
accuracy parameter %(n; Ep(n)). Similarly, we can con- 
sider how much information about p is transferred into 
Pq. We characterize this by considering how much infor- 
mation about n ■ & of p can be obtained by the measure- 
ment of n ■ a on /5q. We can regard the measurement of 
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n ■ a on pq as the measurement described by the POVM 
Eq(«) = {E k M' k f H+;n)M' k , E k Mj}P(-,n)M' k } on 
p. We thus characterize the amount of information which 
is transferred from P to Q by the accuracy parame- 
ter x( n ; Eq(ti)). For mathematical convenience, we use 
ep(n) = e(n;Ep(n)) and £q(»i) = e(n;EQ(n)), instead 
of x( n ; Ep(n)) and x( n ! Eq(ti)), to derive our no-cloning 
inequality. The amount of information about p which re- 
mains in P is characterized by ep(n) averaged over all 
directions, and the amount of information about p which 
is transferred into Q is characterized by sq (n) averaged 
over all directions. 

Definition 4 (cloning parameter). We define the 
cloning parameters Cp and Cq as 



C P = 



jU £p(n) ^' CQ "w Q(n) ^- 



(96) 



Since < £p(n) < oo and < £q>(ri) < oo, the cloning 
parameters satisfy 



< Cp < oo, < C Q < oo. 



(97) 



The cloning parameters depend only on po, p env , and U, 
and characterize the performance of the cloning machine 
{[/, po, /5 onv }. The smaller Cp is, the more information 
about p remains in system P, while the smaller Cq is, the 
more information about p is transferred into system Q 
by the cloning machine. For example, if U is the identity 
operator, then Cp = and Cq = oo hold, which implies 
that all information about p is left in system P. On 
the other hand, if U describes the swapping operation 
between P and Q, then Cp = oo and Cq = hold. 
For intermediate cases between them, the following no- 
cloning inequality between Cp and Cq can be derived 
from trade-off relation (|62|) . 

Theorem 12 (no-cloning inequality): 



CpCq > 



Proof. It can be shown that there exists a POVM 
E(n, n') — {E(i,j;n,n')}, with i,j — ±, satisfying 



((P(i; n) (g> I ® <g> P(j\ n') ® I)Up ®p ® p cnv U f 

= tr [E(iJ;n,n')p 



tr 



(99) 

We can also show that Ep(n) and Eq(n') are its 
marginal POVMs. From inequality (|89|) and the trade-off 
relation we obtain 



ep(n)e Q (n') > e(n; E(n, n'))e(n'; E(n, n')) 

>l-(n-n'f. ( ' 

Averaging (|100|) over all directions and using 



'\n\ = l J\n'\ = l 



£p (n)s Q (n')^^ = CpC Q , (101) 



d 3 n d 3 i 



/-(»■» T)- 47r 



47T 



2 

3' 



(102) 



we obtain 

Inequality (|9~8|) represents the trade-off relation be- 
tween the information remaining in the original system 
P and the information transferred to the reference sys- 
tem Q. The impossibility of achieving Cp = Cq = 
implies the no-cloning theorem. Note that if Cq — > 0, 
then Cp — > oo, which implies that if a cloning machine 
transfers all of the information about p into system Q, 
then no information can be left in system P. 



D. Quantum State Tomography 

We next apply o ur framework to quantum-state tomog- 
raphy [H, [H, HH, H3] • As shown in Sec. VI, charac- 
terization of the measurement accuracy by the accuracy 
matrix is closely related to the asymptotic accuracy of 
the maximum-likelihood estimation which is considered 
to be the standard scheme for quantum-state tomogra- 
phy. 

We first consider the standard strategy to estimate the 
three components of Bloch vector sq- We divide TV iden- 
tically prepared samples into three groups in the ratio 
1:1:1, and measure a x for the first group, a y for the 
second group, and & z for the third group. As N increases, 
this scheme becomes asymptotically described by POVM 
consisting of six operators: 

Ei = \{I + a x ), E 2 = Ul-a x ), 
6 6 

E 3 - \(I + & v ), E A = i(i - a y ), (103) 

E b = hl + a g ), E 6 = \(I-a z ). 
6 6 

We can reconstruct the quantum state by quantum-state 
tomography and hence reconstruct the probability distri- 
butions in all directions. In fact, the accuracy matrix for 
the standard tomography (|103p is given by 



X(E) 




1 



(104) 



which attains the upper bound of the inequality 
Sp(x(E)) < 1. This expression manifestly shows that the 
reconstructive subspace of the standard quantum state 
tomography is R 3 and that the accuracy of the tomogra- 
phy is optimal and symmetric in the sense discussed in 
Sec. HI B. 

We next consider the minimal qubit tomography. 
Rehacek et al. have shown that the following four mea- 
sured probabilities are just enough to estimate the Bloch 
vector [H: 



1 



(1 + Ofc.*) (k = 1,2,3,4), 



(105) 
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where 



ax = 



a 3 




a 2 = 



— 




(106) 



The minimal qubit tomography is also optimal and sym- 
metric, in the sense that the corresponding accuracy ma- 
trix is again given by [[TM]) . Note that the POVM E 
satisfying V(E) = K 3 can be regarded as tomographi- 
cally complete (38j . 



VI. MAXIMUM-LIKELIHOOD ESTIMATION 
AND THE FISHER INFORMATION 

In this section, we point out a close connection between 
the accuracy matrix and the Fisher information [34], |35| . 
We consider the quantum measurements described by 
the POVM E = {E k } for each of N (< oo) sam- 
ples prepared in the same unknown state p. Note that 
Ek = Tk(j + Vf. • <x). Our task is to estimate the Bloch 
vector So by maximum-likelihood estimation. We denote 
s* as the maximum-likelihood estimator of Sq from N 
measurement outcomes. 

The asymptotic accuracy of maximum-likelihood esti- 
mation is characterized by the Fisher information. In our 
situation, the Fisher information takes the matrix form 
given by 



E 



d 2 \nf k (s) 



d(s) i d(s) j 



k qk 



or equivalcntly, 



(107) 
(108) 



Note that I is a 3 x 3 positive and Hermitian matrix, and 
that the support of I coincides with that of %(E). 

Focusing on a particular direction n, we can reduce 
the Fisher information content to 



I(n) 



1 



(109) 



The greater the Fisher information, the more informa- 
tion we can extract from the measurement outcome. In 
the case of I(n) = 0, the variance of the estimator n ■ s* 
diverges, so we cannot gain any information about the 
probability distribution in direction n. This is the case 
of n not being in any reconstructive direction. 

Replacing qk by in the Fisher information (|10T[) or 
(|108|) . we can obtain the accuracy matrix in Eq. (fT2|) or 
(frj|. Note that is the average of q k over the entire 
Bloch sphere. The trade-off relations ([52]). JEH, and (J7SJ) 



can thus be interpreted as the trade-off relations between 
the asymptotic accuracy of the maximum-likelihood es- 
timation of the probability distributions of observables. 
A finite number of samples only gives us imperfect in- 
formation about the probability distribution of an ob- 
servable for an unknown state. As we have shown [23[, 
this imperfection further deteriorates in the case of si- 
multaneous estimation due to the noncommutability of 
the observables. 

Figure 2 shows the results of simulations for the value 
of the maximum-likelihood estimators p(-\-;n A )* (red 
curves) and p(+; ng)* (blue curves) in the the case of an 
optimal nonideal joint POVM which satisfies the equality 
in ((62]) or ([63]) with 



Xa = 1/10 = 0.10, Xb = 36/37 ~ 0.97 



and 




n B 




I + a x 



(110) 



(HI) 




500 1000 
N 



1500 2000 



FIG. 2: (Color online) Maximum-likelihood estimators 
p(+;tia)* (red curves) and p(+;ns)* (blue curves) for the 
case of 6 = 7r/6, \a = 0.10, and \b — 0.97. The abscissa 
indicates sample number N and the ordinate indicates the 
value of the estimators. The number of simulations is 20. 

Let us next consider a simple estimation scheme by di- 
viding N prepared samples into two groups in the ratio 
£ : 1 — £(0<£<1) and performing a nonideal measure- 
ment of tia ■ tr by the POVM {E(±; ua)} with accuracy 
Xa for the former group, and similarly we perform a non- 
ideal measurement of tib • <r by the POVM {E(±; Ub)} 
with accuracy xb for the latter group (see also example 
3 in Sec. II). This measurement can asymptotically be 
described by the POVM 

V = {ZE(±n A ),(l-0E(±n B )}, (112) 

whose accuracy matrix is 

X(E) = £xAn A nl + (1 - OxbUbuI. (113) 
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From Eq. ([40|) in Sec. II, the accuracy parameters in di- 
rections and rig are given by 

X (n A ; E) = £xa, x{n B ; E) = (1 - Oxb, (114) 
and thus 

X (n A ;E)+x(ns;E) < 1. (115) 

We can therefore conclude that a simultaneous measure- 
ment has the advantage over this simple method in that 
the former can access the domain \A + Xb > 1 f° r 
7^ 7r/2, i.e., domain Q in Fig.l. 

VII. SUMMARY AND DISCUSSION 

Projection measurements cannot always be imple- 
mented experimentally. This raises the question of how 
accurately we can obtain information about observables 
from a given imperfect measurement scheme. To quan- 
titatively characterize such measurement accuracy, we 
have introduced the 3x3 accuracy matrix x(E), with 
E = {Ek} being the corresponding POVM. 

We have considered the accuracy matrix of the most 
general class of measurements of a qubit system: general- 
ized simultaneous measurements including nonideal joint 
measurements and quantum-state tomography. From the 
outcomes of generalized simultaneous measurements, we 
can obtain information about more than one observable. 

In terms of the accuracy matrix, we have defined ac- 
curacy parameter x(n; E) and error parameter e(n; E) 
for a direction of n corresponding to the observable 
n ■ &. These parameters satisfy < x( n i E) < 1 and 
< e(n;E) < oo. If x( n ; E) = 1, or equivalently 
e(n; E) = 0, the measurement is equivalent to the pro- 
jection measurement of n ■ <x. On the other hand, if 
x(n; E) = 0, or equivalently e(n; E) = oo, we cannot ob- 
tain any information about the measured system by this 
measurement. 

The accuracy matrix and accuracy parameters give us 
information about observables for which we can recon- 
struct the probability distribution from the measured dis- 
tribution {qk}, where qk = tr(pEk). In fact, we can re- 
construct the probability distribution of observable n ■ <x 
if and only if x(n; E) ^ 0, or equivalently e(n; E) < oo. 
In other words, the direction n is a reconstructive direc- 
tion if and only if n e V(E), where the subspace V(E) of 
R 3 is spanned by the eigenvectors of x(E) corresponding 
to nonzero eigenvalues. 

The main results of this paper are trade-off relations 
(l62j) . (|63l) . and (|76|) between the accuracy parameters 
and the error parameters. We can interpret them as the 
uncertainty relations between measurement errors in gen- 
eralized simultaneous measurements; the more informa- 
tion we obtain about an observable, the less information 
we can access about the other noncommuting observable. 
Trade-off relation (|6"2")1 can also be interpreted as the un- 
certainty relation between the measurement error and 



back-action of measurement as formulated in inequal- 
ity (USD. 

The new no-cloning inequality in (|98[) is derived from 
the trade-off relations. To derive this, we have introduced 
the cloning parameters Cp and Cq , where P indicates the 
system to be cloned and Q indicates the blank reference 
system. Let p be the pre-cloned state of system P. After 
a cloning operation, all the information about p remains 
in system P if and only if Cp = 0, and the information 
about p is completely transferred to system Q if and only 
if Cq = 0. The impossibility of attaining Cp = Cq = 
implies the no-cloning theorem. The condition of the 
equality in our no-cloning inequality (|98[) has yet to be 
understood. 

We have also applied the trade-off relations to analyze 
the efficiency of quantum-state tomography. The accu- 
racy matrix of the standard qubit-state tomography or 
the minimal qubit tomography is given by x(E) = Is/3 
with 1% being the 3x3 identity matrix, which implies 
that the efficiency of quantum-state tomography is opti- 
mal and symmetric. 

We have pointed out a close relationship between the 
accuracy matrix and the Fisher information. We have 
also shown that the trade-off relations can be interpreted 
as being those concerning the accuracy of the maximum- 
likelihood estimators of the probability distributions of 
noncommuting observables. 

While we focus on the spin- 1/2 system in the present 
paper, many results can be generalized for higher- 
dimensional systems. We conclude this paper by out- 
lining such generalization. 

In the case of a <i-dimensional system (d > 3), the 
parametrization of the Hcrmitian operator E is given by 

E = r(I+ \fd~Tv ■ A), (116) 

where r is a real number, v is a d 2 — 1-dimcnsional real 
vector, and A = (Ai, A2, • • • , \,p-i) is the elements of the 
Lie algebra of SU(d) satisfying tr(Aj) = and tr(AtAj) = 
dSij with dij being the Kronecker delta. The necessary 

and sufficient condition for E to be a positive operator 
is given by r > and S m (v) > (m = 2, • d) , where 
S m (v) is an mth-degree polynomial for v [42, [43] . The 
condition for m = 2 is given by S^u) = d(d — l)r 2 (l — 
I f 1 2 ) /2 > 0, which is equivalent to 

M<1. (117) 

For m = 3, S^(v) is given by 

S 3 (v) = ^d(d~l)(d-2)r 3 
6 

\ ijk J 

(118) 

where dijk is defined as {A;,Aj} = 2SijI + J2kdijk^k 
with {A, B} = AB + BA. We note that if E is a rank-1 
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projection operator, then \v\ = 1. However, the Hermi- 
tian operator E with \v\ = 1 is not necessarily a positive 
operator. 

The accuracy matrix for a ci-dimensional system as- 
sumes the same form as Eq. (|14[) using parametriza- 
tion (|116p . In this case, x(E) is a d 2 — 1 square ma- 
trix. Moreover, we can define the accuracy parameter 
and the error parameter according to Eqs. (|27| and (|28|) . 
respectively. Using condition (|117[) . we can derive trade- 
off relations ([62]) and (|63|) for a <i-dimensional system. In 
this sense, the trade-off relations serve as universal un- 
certainty relations holding true for all finite-dimensional 
systems. 

However, bounds of trade-off relations (|62|) and ([63]) 
would not necessarily be able to be reached for d > 3, 
because r > and \v\ — 1 are not sufficient for positivity 
of the POVM. Moreover, while the accuracy parameter 
x( n 'i E) for d = 2 characterizes the measurement accu- 
racy of spin observables n ■ <x, the accuracy parameter 
x( n 'i E) for d > 3 cannot characterize the measurement 
accuracy of, for example, the spin-d observable J z ; it only 
characterizes the accuracy of a rank-1 projection opera- 
tor. Therefore the results of this paper based on x( n i E) 
cannot be applied straightforwardly for d > 3. A full 
investigation of this problem is underway. 



APPENDIX A: PROOF OF THEOREM 10 

We prove the case of a = A. For simplicity of no- 
tation, we define that E\ = E(+,+), E 2 = £(+,—), 
i?3 = E(— , +), and £4 = E(— , — ). The accuracy matrix 
is given by 



X (E) = ^>K| 2 , 



(Al) 



k=l 



and the accuracy parameter in direction n A is 

XA=(n A -( X (E)- 1 )n A y 1 . (A2) 
The marginal POVM E^ is 

E A (+) = (n+r 2 )(i+ nVl+ / 2V2 - ( r), (A3) 
V ri + r 2 I 



tr ! \ 1 , ^ I t , r3Vs + r4V4 
Ea{-) = (r 3 +r 4 ) ■ er , 

and the marginal accuracy matrix is 



x(Ea) = x\^ A n 



where 



X A 



Ai 



\rw x + r 2 v 2 \ 2 , |r 3 i;3 + r 4 v 4 \ 2 



ri + r 2 



T3 + ri 



(A4) 



(A5) 



(A6) 



Our objective is to show that xa = x' A - For simplicity, 
we introduce the notation 




(A7) 



We can then write x(E) as 



-(E) = a k a k = (Xx, Xy, Xz 



(A8) 



fe=l 



where 



Xx = ^2 x k a k , Xy = ^2 Vka k , Xz = ^ z k a k . (A9) 



k=l 



k=l 



fe=l 



Using Eq. (|A9[) . we can calculate the determinant of 
X(E): 



detx(E) = ^x k a k J ■ ^y k a k x ^z k a k 

\k=l J \k=l k=l / 

1 2 



Y [ a k ■ (ai x am)]' 



k<l<m 



(A10) 



On the other hand, the cofactor matrix of x(E) is 

/ (Xy X Xz)x (Xy X Xz)y (Xy X Xz) z 
X(E) = (Xz X Xx)x (Xz X Xx)y (Xz X Xx) z 
V (Xx X Xy)x (Xx X Xy)y (Xx X X y ) 2 

Therefore 

x(E),j = ^2(a k x at)i(a k x a ; )j. 

The inverse matrix is given by 

, f^-is Efc<;( Q fc x ai)i(a k x ai)j 

(X(E) - — -2 ■ 

T,k<l<m K ■ ( a i x a m )] 



(All) 
(A12) 

(A13) 



Noting that 
(nvi + r 2 v 2 ) ■ x(E)(nvi + r 2 v 2 ) 

= ^2^2( r i v i + r 2'"2)i(afc x ai)i(a k x aj^riDi + r 2 v 2 ) 

= [(n«i + r 2 t; 2 ) • («i x a 3 )] 2 

fc<z 

= (n. + r 2 ) ( [ai • (a 2 x a 3 )] 2 + [a x ■ (a 2 x a 4 )] ' 



and 



n A 



r 1 v 1 + r 2 t> 2 
IriUi + r 2 t> 2 | 



(A14) 



(A15) 



14 



we obtain 



n+r 2 [ai ■ (02 x a 3 )] 2 + [ai • (a 2 x a 4 )] 2 
|ri«i + r 2 f 2 | 2 detx(E) 

(A16) 

Similarly, we can show that 

r 3 + r A [a 3 ■ (a 4 x ai)] 2 + [a 3 • (a 4 x a 2 )] 2 



r 3 t> 3 + r 4 v 4 



detx(E) 



Let us dehne 

[ai ■ (02 x a 3 )] 2 + [a a • (a 2 x a 4 ) 



detx(E) 



(A17) 



(A18) 



We can thus conclude 

n A -( X (E)- 1 )n A = t(l-t)ri. 

Therefore 

1 



XA 



n A ■ ( X (E)-i)n A 
1 

V(n + r 2 ){r 3 + r 4 ) 
1/1 1 



|ri«i + r 2 t? 2 | 2 , \r 3 v 3 + r i v i \ 2 



r\ + r 2 



r 3 + r4 



which is our objective. 



(A22) 



(A23) 



_ [03 ■ (04 x ai)] + [g 3 ■ (g 4 x g 2 )] 

det X (E) ' lAi9j 

i = n + r 2 = 1 - (r 3 + r 4 ), and l/?7 = |n«i + r 2 v 2 | 2 = 
|r 3 u 3 + r 4 w 4 | 2 . Noting that A + Y = 1 and 

n A • (xCE)- 1 )^ = Xtrj = Y(l - tjjj, (A20) 
we obtain 

X + t = l. (A21) 
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